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Abstract 

We analyze a number of "nearly exactly" solvable models of electronic spec- 
trum of two-dimensional systems with well-developed fluctuations of short 
range order of "dielectric" (e.g. antiferromagnetic) or "superconducting" 
type, which lead to the formation of anisotropic pseudogap state on certain 
parts of the Fermi surface. We formulate a recurrence procedure to calculate 
one-electron Green's function which takes into account all Feynman diagrams 
in perturbation series and is based upon the approximate Ansatz for higher- 
order terms in this series. Detailed results for spectral densities and density 
of states are presented. We also discuss some important points concerning 
the justification of our Ansatz for higher-order contributions. 

PACS numbers: 74.20.Mn, 74.72.-h, 74.25.-q, 74.25.Jb 
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I. INTRODUCTION. 



Recently there was an upsurge of interest in observation of the pseudogap in the spec- 
trum of elementary excitations of underdoped high-temperature superconductors (HTSC) 
|l]||. These anomalies were observed in a number of experiments, such as the measurements 
of optical conductivity, NMR, inelastic neutron scattering and angle-resolved photoemis- 
sion (ARPES) (cf. review in Q). Probably most striking evidence for the existence of this 
unusual state were obtained in ARPES-experiments which demonstrated essentially 

anisotropic changes in the spectral density of current carriers in rather wide temperature 
region in normal (nonsuperconducting) phase of these systems (cf. review in 0). A re- 
markable anomaly observed in these experiments was that the most significant changes (in 
comparison with the usual Fermi-liquid behavior) in spectral density were seen close to the 
point (71", 0) in the Brillouin zone, while there were no changes at all in the direction of zone 
diagonal (point (tc,tt)), which in fact demonstrates the "destruction" of the Fermi surface 
around the point (tt, 0) and conservation of Fermi-liquid behavior in the direction of zone 
diagonal. In this sense it is usually claimed that the pseudogap symmetry is of the "rf-wave" 
type, the same as the symmetry of superconducting energy gap in these compounds 
At the same time the mere fact that these anomalies exist at the temperatures much higher 
than superconducting transition temperature, as well as in the underdoped (non optimal) 
region of carrier concentrations, may signify some other nature of these anomalies, not con- 
nected directly with Cooper pairing. Typical phase diagram of HTSC-system is shown in 
Fig.l. 

There are rather many theoretical papers now attempting to interpret the observed 
anomalies. We may classify these attempts in two main schools of thought (scenarios). 
One is based upon the idea of Cooper pairs (fluctuation) formation at temperatures higher 
than the usual superconducting transition temperature Another assumes that the 

pseudogap phenomena are determined mainly by fluctuations of antiferromagnetic (AFM) 
short range order [|]-[12 ■ 

Rather long ago one of the authors of the present paper had proposed an exactly solvable 
model of the pseudogap formation in one-dimensional system due to well developed fluctu- 
ations of short range order of charge density (CDW) or spin density (SDW) wave fi"3|-|T7 . 



Recently this model has attracted some interest in connection with attempts to understand 



the pseudogap state in HTSC-cuprates [11,12,18-20]. In particular in Refs. [11,12 an impor- 



tant generalization of this model was formulated for the case of two-dimensional electronic 
system with random field of spin fluctuations (AFM short range order). In this model of 
"hot spots" on the Fermi surface |p!IU2| , using the formal scheme of Refs. [|l5HT7l|, they 



obtained rather detailed description of pseudogap anomalies for the case of large enough 
temperatures ("weak pseudogap" region in Fig.l). In Refs. |T9|j20|] a simplified variant of 
this model |T3]fn}j , corresponding to the limit of very large correlation lengths of fluctuations 
of short range order, was used to describe the pseudogap state determined by well devel- 
oped fluctuations of superconducting (SC) short range order. In a recent paper ET1] this 
(over) simplified model was used to derive Ginzburg-Landau expansion (for different types 
of Cooper pairing) for the system with strong fluctuations of CDW(SDW,AFM)-type, in 



the framework of "hot patches" model proposed in this work. At the same time, in Ref. [22 



dedicated to rather detailed review of the model proposed in Refs. [|I"3|-|T7H, a certain error 
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was noted in early papers [1~5|-|I7| in the analysis of the case of finite correlation lengths of 
fluctuations of short range order. In Ref. [12] it was claimed that this error is rather insignif- 
icant, especially in two- dimensional model of "hot spots", which is of the main interest for 
the physics of HTSC-cuprates. 

The aim of the present paper is to perform an analysis of number of important aspects 
of this "nearly exactly" solvable model, mainly for two-dimensional case. We shall consider 
both the case of short range order fluctuations of CDW(SDW,AFM)-type in the "hot spots 
" model fTl,12|, as well as the possible applications of this model in the scenario of supercon- 
ducting fluctuations [[7|, |iT?pU| (SC short range order), in particular for the most interesting 
case of (i-wave pairing. Besides the general discussion of reliability of the formal scheme 
of Refs. |T3 Hr7ip!Il , p!2l , we shall perform detailed calculations of the spectral densities and 
one-electron density of states both for the "hot spots" model [0],O], and in the scenario of 
SC-fluctuations. 



II. "HOT SPOTS" MODEL. 

1. Model description and "nearly exact" solution for the Green's function. 



The model of "nearly-antiferromagnetic" Fermi-liquid p3| , |24|| is based upon the picture 
of well developed fluctuations of AFM short range order in a wide region of phase diagram 
shown in Fig.l. In this model the effective interaction of electrons with spin fluctuations is 
described via dynamic spin susceptibility % q (u;), which is determined mainly from the fit to 
NMR experiments [ p4|j : 

K//(q, w) = g\M « 1 , C ^/^L (1) 



i+e(q-Q) 2 -*- 



'sf 



where g is coupling constant, ^-correlation length of spin fluctuations, Q = (n/a,n/a)- 
vector of antiferromagnetic ordering in insulating phase, w s y- characteristic frequency of 
spin fluctuations, a-lattice spacing. 

As dynamic spin susceptibility Xqi 00 ) nas peaks at the wave vectors around (ir/a, ir/a) 
there appear "two types" of quasiparticles — "hot quasiparticles" with momenta in the 
vicinity of "hot spots" on the Fermi surface (Fig. 2) and "cold" quasiparticles with momenta 
on the other parts of the Fermi surface, e.g. around diagonals of the Brillouin zone \p x \ = \p y \ 
TT|,|T2"| • These terms are connected with the fact that quasiparticles from the vicinity of "hot 



spots" are strongly scattered on the vector of the order of Q by spin fluctuations (|l|), while 
for quasiparticles with momenta far from "hot spots" this interaction is relatively weak. 
In the following we shall consider the case of high enough temperatures when nT ^> uj s f 



which corresponds to the region of "weak pseudogap" in Fig.l [p] , p^| . In this case spin 



dynamics is irrelevant and we can limit ourselves to static approximation: 

^M = A, i + { .£-0). < 2 > 

where A-effective energy parameter, which in the model of AFM fluctuations can be written 
as IOI: 
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A 2 = g 2 Tj2x^ m )=9 2 <S*>/3 (3) 

mq 

where Si-spin on the lattice site {Cu ion in highly conducting Cu0 2 -pl3Jie for HTSC- 
cuprates). In the following we shall consider A (as well as £) as some phenomenological 
parameter, determining the effective width of the pseudogap. 

We can greatly simplify all calculations if instead of (y) we use another form of model 
interaction (cf. analogous model in ||): 

K//(q) = A V +£-<?.)» {->+£-<?,)» (4) 

where A 2 = A 2 /4. In fact (f|) is qualitatively similar to (0) and differs from it very slightly 
in most important region of |q — Q| < 

Consider the first order in V e ff correction to electron self-energy shown in Fig. 3: 

£( gn p) = ^> e// (q) * (5) 

q ten Sp+q 

Main contribution to the sum over q comes from the region close to Q = (ir/a, it /a). Then 
we can write: 

£p+q = £p+Q+k ~ £p+Q + V p+Qk (6) 
where fp + Q = 8 |^ +Q , and performing integration over k, we get: [] 

A 2 

ie n - 4 p+ q + U^p+qI + \ v p+Q\) KS Wn £ n 

where k = 

The spectrum of "bare" (free) quasiparticles can be taken in the form [0.0 : 

£ p = — 2t(cosp x a + cospyd) — At cosp x acosp y a (8) 

where t-nearest neighbor transfer integral, t'-second nearest neighbor transfer integral on 
the square lattice, [i is chemical potential. In the analysis of real HTSC-systems in Refs. 



n,|12|| it was assumed e.g. for YBa 2 Cu 3 6+ s that t = 0.25eV, t = — 0.45t , while fi was 



fixed by hole concentration. Below we shall see that it is of considerable interest to study 
our problem for different relations between t and t . 



x In Refs. [ll,nj another but similar in spirit to (Q) form of effective interaction was used: 

' rCk-f r a +ki ' where k \\(±) 

that result analogous to (0) takes the form: 



V e //(k) = A 2 2 | 2 jMrcz i where k\uj_\ -projection of k parallel (perpendicular) to v p+ q, so 



is - £ p+ q + i\v p+ Q\Ksigne n 
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Consider second order corrections to self-energy shown in Fig. 4. Using (^J) we obtain: 

dki f <ik 2 k k k k 



no) = A 



-2 / T 2 „2 I h.2 K 2 _|_ h.2 K 2 i h.2 K 2 i p 



l l 



ie n - £ p+Q - fp+Q^Lr - ^p+Q^lj/ ~ £p ~ ^(^la + &2:r) ~ fp^l?, + fc 2j ,) 

1 



W ~~ / —2 / -7r2 ^.2 _|_ p ^.2 _|_ 1-2 ^.2 _i_ p .,2 _|_ JU2 

«/ /I »/ /i /u la; lj/ 2a? 2?/ 

l l 

- £ P+Q - fp+Q^ix - ^p+Q^iy ie n - £ P - v*{k lx + k 2x ) - v%(k ly + k 2y ) 

1 



(9) 



ten - £p+Q - V* +Q k 2x - V V p+Q k 2y 



(10) 



where we have used the spectrum (§), from which, in particular, it follows that £ P +2Q = £ P , 
v p +2Q = v p for Q = (n/a, ir/a). If signs of and ^p+Q, and of and v^ + q are the same, 
we can see that integrals in (|9]) and (|10D are fully determined by poles of Lorentzians, which 
determine the interaction with fluctuations. After elementary contour integration we obtain: 



1 1 

s(a) = m = Wn - e P+Q + i(K^\ + K +Q \h? ten - e P + Wh\ + \W* 

Here and in the following we assume e n > 0. It is easy to convince ourselves that in 
case of velocity projections of the same sign we can similarly calculate contributions of any 
higher order diagrams. Accordingly, the contribution of an arbitrary diagram for electron 
self-energy of N-th order in interaction (0) has the form: 

2N-1 1 

s^M = A 2Ar II - — rr (12) 

where £j = £ p +q and Vj = I^p+qI + I u p+q| f° r °dd j and £j = £ p and Vj = \Vp\ + \v*\ for even 
j. Here rij is the number of interaction lines, enveloping j-th Green's function in a given 
diagram. 

In this case any diagram with intersecting interaction lines is actually equal to some 
diagram of the same order with noncrossing interaction lines. Thus in fact we can consider 



2t 



In the model of V e ff used in Ref. [11,1^], for the case of v p v p+ Q > 0, in a similar way we get: 

E(o Eb A [j £n _ ^ p+Q + i|v p+Q | K ] 2 ie n - ^p+Q + i2|v p |(| cos(/)| + | sin^|)K 
where i^-angle between v p and v p+ q. 
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only diagrams with nonintersecting interaction lines, taking into account diagrams with 
intersecting lines introducing additional combinatorial factors into interaction vertices. This 
method was first introduced (for another problem) in a paper by Elyutin |25j and used for 
one-dimensional model of the pseudogap state in Refs. |T5Hl7| • 

As a result we obtain the following recursion relation for one-electron Green's function 
(continuous fraction representation) flTEHTi 



Go 1 ' 



£n£ P ) — £i(e n £p) 



(13) 



S,(^p) = A 2 - - f - .J {k) = , p . (14) 

where £ fc = £ P+Q and v k = \Vp +Q \ + I^p+qI for odd values of k, £ fc = £ p and v k = \v*\ + 
for even k. Combinatorial factor: 



v(k) = k (15) 



corresponds to our case of commensurate fluctuations with Q = (7r/a, 7r/a) [15|]. It is not 
difficult to analyze also the case of incommensurate fluctuations when Q is not "locked" 
to the period of inverse lattice. In this case diagrams with interaction lines enveloping odd 
number of vertices are significantly smaller than diagrams with interaction lines enveloping 
even number of vertices. Thus we must consider only these last diagrams ||13|-0| . Recurrence 



relation has the same form fljj), but diagram combinatorics and factors v(k) change [15| : 



v(k) = { f 5° r ° dd \ (16) 
I I tor even k 



In Refs. |TT| , p!2|] a spin-structure of effective interaction within the model of "nearly an- 



tiferromagnetic" Fermi-liquid was taken into account (spin-fermion model of Ref. [0). 
This leads to more complicated combinatorics of diagrams in the commensurate case with 
Q = (it /a,ir/a). Spin-conserving part of the interaction gives formally commensurate com- 
binatorics, while spin-flip scattering is described by diagrams with combinatorics of incom- 
mensurate type ( "charged" random field in terms of Ref. [T2|] ) . As a result the recurrence 
relation for the Green's function is again of the form of(0), but the combinatorial factor 



v{k) is now [ ll|,|T2[| : 



v t k \ = { ^f- for odd k n 7) 

1 I for even k 

As we noted above, the solution of the form of ( |14|) can be obtained only in the case 
of coinciding signs of velocity projections Vp + Q(fp + q) and fp(fK). Below we shall analyze 
situations when this is really so. In case of different signs of these projections integrals of the 
type of (0) and fllUD , corresponding to higher order corrections, can not be calculated in such 
a simple way as above, because contributions from the poles of the Green's functions become 
relevant. In this case, instead of simple answers like QTTJ) rather cumbersome expressions 
appear and, even more importantly, disappears the fundamental (for our method) fact of 
equality of wide range of diagrams with crossing and noncrossing interaction lines, which 
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actually allows us to classify higher order contributions and obtain an "exact" solution 
(14). This problem is important only for the case of finite correlation lengths of fluctuations 
£ = while in the limit of £ — - > oo(k — > 0) the exact solution for the Green's function 
is independent of velocities v p and v p+ q and can be easily obtained in analytic form by 
methods of Refs. |pi|,|l4"| (cf. also [I2]|). In one-dimensional model considered in Refs. JD^D]] 
the signs of corresponding velocity projections are always different (these correspond to 
"left" and "right" moving electrons). This problem was stressed in a recent paper P2"| . In 
Appendix A we analyze these difficulties in detail for one-dimensional case and show that 
the "Ansatz" of the type of ( |P2| ) used in Refs. |T5|-[T7| for the contributions of higher order 
diagrams as well as solution of the form of ([H]) in fact give us very good approximation to 
an exact solution even in the case of velocity projections of different signs. Obviously this 
solution is exact in the limits of £ — > oo(k — > 0) and £ — > 0(k — > oo), and guarantees rather 
good (qualitatively) description in the region of finite correlation lengths. 



2. Analysis of the "bare" energy-spectrum. 

For the "bare" energy spectrum (§) we can easily find conditions (relations between t,t' 
and fi), when solution (IT1) is exact. First of all, let us define the region of parameters t,t' 
and n, when "hot spots" on the Fermi surface exist, i.e. the condition of existence of points 
connected with vector Q = (n/a, ir/a). If p = (p x ,p y ) is the position of a "hot spot" on the 
Fermi surface, then p + Q = (p x + ir / a, p y + ir / a) must also belong to the Fermi surface, so 
that for the spectrum (H) we must have: 

— 2t(cos p x a + cos p y a) — At' cos p x a cos p y a — \i = 

2t(cosp x a + cospyo) — At' cosp x a cosp y a — /i = (18) 



Then the condition of existence of "hot spots" becomes: 

2 n„a = 

At 



cospyd = —cosp x a and cos 2 p x a = — (19) 



Thus, the "hot spots" on the Fermi surface exist if: 



< < 1 (20) 
~ At' ~ K ' 

Define now the region of parameters t,t' and /i where (0) is exact requiring the positivity 
of products fpfp + Q and fpfp+Q- We have: 

= — — = It sin p x a + At' sin p x a cos p v a 
P dp x 



v. 



P dp y 



2t sin p y a + At' sin p y a cos p x a 



v p v p+q = IQt' 2 sin 2 p x a[cos 2 p y a - (^) 2 ] 

v v p v y p+Q = 16t' 2 sin 2 p y a[cos 2 p x a - (^) 2 ] (21) 



S 



It is easily seen that for the existence of points on the Fermi surface, where the projections of 
velocities are of the same signs, it is necessary to fulfil \t'/t\ > 1/2. We are mainly interested 
in the region around the "hot spots" , where with an account of(pT9|) we have: 



« +Q = « +Q = 4t 2 (l-| 7 )(^-l) (22) 
Thus the projections of velocities in "hot spots" have the same sign if: 



fit' 
7 



> 1 (23) 



The same condition obviously guarantees the validity of inequality v p v p+ q > which is 
necessary for exactness of (|T4] ) in the model of Refs. flTT|J12"|| . 

In Fig. 5 we show (dashed region) the region of parameters where "hot spots" exist 
(0 < /i/4t' < 1) as well as the region where velocity projections in their vicinity are of the 
same sign (fit' > 1). In Fig.6 we show the Fermi surfaces, defined by the "bare" spectrum 
for different values of chemical potential fi (band fillings) when these conditions are 
either satisfied or not. 



3. Spectral density and density of states. 
Consider one-electron spectral density: 

A(Ep) = —ImG R {Ep) (24) 

7T 

where G R (Ep) is retarded Green's function, obtained by the usual analytic continuation of 
flTjl) to the real axis of energy E. In Fig. 7 we show energy dependencies of A(Ep) obtained 
from fll3D,(|i~4]) for different variants of combinatorial factors (|i~5[ ) , ((T6| ) , (|I~7 D . For t'/t = —0.6 
and fi/t = —1.8 < t/t' = 1.666 projections of velocities in "hot spots" are of the same sign 
and relation ( |14"|) defines Green's function exactly. We can see that in incommensurate case 
(|I~6|) (Fig.7(a)) as well as for combinatorics of spin-fermion model ( |TTD (Fig. 7(c)) the spectral 
density at the "hot spot" demonstrates clearly non Fermi-liquid behavior (for large enough 
values of correlation length £). Note also that for both cases the numerical values of spectral 
density are very close. In the case of commensurate combinatorics (|15|) (Fig. 7(b)) precisely 
at the "hot spot" the spectral density has a single peak and, in this sense, is similar to that 
of the Fermi-liquid even for large values of £. However, even in the nearest neighborhood 
of "hot spot" this spectral density acquires two peak structure ("shadow" band) of non 
Fermi-liquid type for large enough values of £ (see insert in Fig. 7(b)). 

Far from "hot spots" velocity projections are, in general, of different signs even if con- 
dition Q2"Bp is satisfied. Accordingly, the recurrence relation ([TJ]) for the Green's function 
is non exact here. [] However, from discussion in Appendix A it becomes clear that our 



3 At the same time, with the growth of £ more narrow vicinity of the "hot spot" becomes important 
and our approximation is more and more accurate. 
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"Ansatz" fll2|) and solution ( fT4] ) in fact somehow overestimate the role of finite correlation 
length £. There we also provide slightly different variant of solution ( |A11|) , which somehow 
underestimates this role. Inserts in Fig. 7 we show energy dependencies of spectral density 
far from the "hot spot" for different combinatorics of diagrams (|15|), (0), (|0)- 

For comparison we also show the spectral density (for incommensurate case) obtained 
with an "Ansatz of alternating k" ( |A1 1|) which happens to be very close to that obtained 
from flTB|), (|H|). This confirms rather high accuracy of ( |T4D for arbitrary momenta close to 
the Fermi surface. 

In Fig. 8 we show energy dependencies of spectral density for different combinatorics 
(f|),(|D,(|7|) at the "hot spot" in the case of t'/t = -0.4, which, according to Ref. UnjjEZ 



corresponds to YBa2CusO^ + s. For this value of t/t! even at the "hot spot" velocity projec- 
tions are of different signs. However, spectral density (for incommensurate case) obtained 
from the solution with "alternating" k (dashed line in Fig. 8 (a)) is seen to be very close to 
that obtained from (0). This shows that "Ansatz" flT2"| ) and solution (|TJj) produce results 
which are quantitatively close to an exact solution. Let us stress once again that our solution 



(|I~4|) is exact both for £ — > oo and £ — ► 0, while in the region of finite £ it provides rather 
good interpolation. 

Consider now one-electron density of states: 

N(E) = £ A(E, P ) = -- E ImG R (Ep) (25) 
p n p 

which is determined by the integral of spectral density A(Ep) over the Brillouin zone. We 
have seen above that though for some topologies of the "bare" Fermi surface (band fillings) 
we can guarantee the same signs of velocity projections close to the "hot spots" , this is not 



so in general case far from "hot spots", so that solution (|14f) , based upon our "Ansatz" ( jl2| ) 



is only approximate. Accordingly, the use of (|14[) to calculate the density of states with ( |25|) 



leads also to a kind of approximation. In Fig. 9 we show densities of states obtained from 
(|Ilf),(|l4|),([2"5|), with the use of the "bare" spectrum (|8|) for different diagram combinatorics 



©,(0),©, for the values of t'/t = -0.4 (Fig.9(a)) and t'/t = -0.6 (Fig.9(b)). We can see 
that for t'/t = —0.4 there appears some dip in the density of states (pseudogap). This dip is 
relatively weakly dependent on the value of correlation length £ (see insert in Fig.9(b)). If 
the band filling is appropriate and the Fermi level [i is somewhere in this energy region there 
are also "hot spots" at the Fermi surface. For t'/t = —0.6 the region of existence of "hot 
spots" is rather wide, but the pseudogap in the density of states is practically unobservable. 
We can see only the obvious smearing of Van-Hove singularity which is present in the absence 
of fluctuations. 



III. MODEL OF "SUPERCONDUCTING" FLUCTUATIONS. 



4- Model description and solution for the Green's function. 

Pseudogap phenomena can also be probably explained using the ideas of fluctuation 
Cooper pairing at temperatures higher than superconducting transition temperature T c 
us cons ider the simplest possible model approach to this problem. In Fig. 10 (a) 
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we show the self-energy diagram of the first order in fluctuation propagator of Cooper pairs 
for T > T c . Anticipating the possibility of both usual s-wave and d-w&ve pairing, charac- 
teristic of HTSC-systems, we introduce the pairing interaction of the simplest (separable) 
form: 

V(p,p') = -Ve^e^') (26) 

where <fi is polar angle determining the direction of electronic momentum p in the plane, 
while for e(0) we assume model dependence pB]f27|] : 



= J 1 s-wave pairing 
| v^2cos(20) d-wave pairing 

Interaction constant V is as usual assumed to be non zero in some energy layer around the 
Fermi surface. Then the self-energy corresponding to Fig. 10 (a) takes the form: 

S(e n p) = V ef f(iuj m c\)G(iuj m - ie n , -p + q) (28) 

mp 

where effective interaction with SC-fluctuations can be written as: 

V eff (iu m q) = - i _ VT G Q (iE n p)G Q {f! m - ie nt -p + q)e 2 (0) (29) 



In the following we assume SC-fluctuations static fpo that in (|33|) we can limit ourselves 



only to the term with uj m = 0. Then effective interaction can be written as: 



where 



: / ; ^o^o.is^ (3i) 



is the usual coherence length of superconductor, A 2 = N ^ F ^2 (N(E F )-density of states 

at the Fermi level Ep). Surely, within the simplest BCS-like model used here, we have 
A w 2n 2 T c (T c / E F ) ~ A (A /_E^) A (where A is superconducting energy gap at 
T = 0) and an obvious problem to explain the experimentally observable scale of pseudogap 
anomalies appears. However, in the following we again consider both £ and A as some 
phenomenological parameters to be determined from experiment on HTSC-systems and not 
from naive BCS-like model. 



4 Static approximation here is valid for ttT 3> u> sc = 8(T — T c )/ir, which is formally analogous to 
the condition of irT S> w s t used in the "hot spot" model above. Here it is well satisfied if the 
system is close enough to superconducting transition 
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Analogously to transition from (|2]) to (f|) we introduce instead of (^) the model inter- 
action of the form: 

V eff ( q ) = -AV(0) J; r \ (32) 



where A = A /4. Quantitatively this is close enough to (130) and leads to great simplifi- 



cation of calculation allowing us to classify contributions of higher order diagrams. In this 
case the first order contribution of diagram in Fig. 10(a) is: 

%£n + £ P + i{\v x \ + \v y \)Ksigne n 
where v x = VfCOS(p,v y = Vpsm(j),K = Second order contribution from Fig.lO(b) is: 

2 f dqi x k f dq\y k f dq\ x k f dq\ y k 



S (2) M = (A 2 e 2 W) : 



7r k -f- qf x J tt k z + qf y J 7i n z + qi x J f qi y 

1 1 
(i£ n + £p~ viqi) 2 ie n - ^ p - v 2 qi - v 2 q 2 



(34) 



where vi = — V2 = vp. In fact we can easily see that in this problem we have practically 
the same rules of diagram technique as in the "hot spot" model, but with combinatorics 
of incommensurate case. This last fact is obvious from the topology of interaction line 
(fluctuation propagator of Cooper pairs) in diagram of Fig. 10 (a) — it is seen that in higher 
orders only those diagrams exist in which interaction lines envelop only even number of 
interaction vertices. The expression of (0) is quite analogous to that of (0), but the signs of 
velocity projections in denominators of Green' functions here are always different: vi = — V2. 
Thus in diagrams of higher orders there appear contributions not only from Lorentzians of 
interaction, but also from the Green's functions. However (in view of discussion in Appendix 
A) we can estimate contributions of higher order diagrams using the "Ansatz" of the type 



of (0), i.e. calculate all integrals e.g. in fl3~4]) , as if velocity projections vi and V2 are of the 
same sign, but in answer just put vi = — V2 = v^. Then we again obtain the recurrence 
relation for the Green's function of the type of (|i~4"|): 

k np ie n - (-l) fc £ P + ikv F K(\ cos 0| + | sm<p\) - ^ k+1 (e n ^ p ) > " 

where v (k) is defined in ([H]). Surely, this relation ( j35|) is not exact, but again it gives exact 
results for the limits of k — > 0(£ — ► oo) and k — > oo(£ — ► 0) and provides rather good 
(quantitatively) interpolation between these limits for the case of finite correlation lengths. 



5. Spectral density and density of states. 

In Fig. 11 (a) we show energy dependencies of the spectral density A(Ep) for one-particle 
Green's function (|24"D, calculated from ( |3"5| ) for different values of polar angle 0, determining 
the direction of electronic momentum in the plane (we take here |p| = pf), for the case of 
fluctuations of <i-wave pairing. It is clearly seen that in the vicinity of the point (7r/a, 0) 
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in Brillouin zone this spectral density is non Fermi-liquid like (pseudogap). As vector p 
rotates in the direction of the zone diagonal the two peak structure gradually disappears 
and spectral density transforms to the typical Fermi-liquid like with a single peak, which is 
narrows as approaches 7r/4. Analogous transformation of the spectral density takes place 
as correlation length £ becomes smaller. 

In Fig. 11(b) we also show the evolution of the product f(E)A(Ep) (where f(E) is fermi 
distribution) which is essentially the parameter measured in ARPES experiments D. Note 
that curves in Fig. 11(b) are quite similar to those obtained in Refs. [|iT||L2|| within the "hot 
spots" model. This picture of Fermi-surface destruction following from this calculations is 
qualitatively shown in Fig. 12 and is very similar to experimental data obtained e.g. in Ref. 
H for Bi 2 Sr 2 CaCu 2 8+ 5. 

In case of fluctuation pairing of s-wave type the pseudogap appears isotropically on the 
whole Fermi-surface and spectral density is non Fermi-liquid everywhere for the case of large 
enough correlation lengths £ of SC-fluctuations. 

In Fig. 13 we present the results of calculations of one-electron density of states with the 
help of ( |35| ) both for the case of s-wave pairing (Fig.l3(a)) and <i-wave pairing (Fig. 13(b)), 
for different values of correlation length of SC-fluctuations. We see that in the case of ci-wave 
pairing the pseudogap is the density of states is naturally not so deep as in s-wave case, 
even for large enough correlation lengths. At the same time it is seen that the pseudogap in 
the density of states in the model of SC-fluctuations is nevertheless much more expressive 
than in the model of "hot spots" due to AFM-fTuctuations. 



IV. CONCLUSION. 

In this paper we have studied "nearly exactly" solvable models of the pseudogap state 
of the electronic spectrum of two-dimensional systems, based upon alternative scenarios of 
its origin — the picture of fluctuations of "dielectric" (AFM, SDW, CDW) type, leading to 
the model of of "hot spots", and the picture of fluctuational Cooper pairing above T c (SC- 
fluctuations). The term "nearly exactly" solvable means that in this approach we can sum 
all Feynman diagram series for one-particle Green's function (and in fact also for two-particle 
Green' function |I6| , [17| ), using for the higher order diagrams an approximate "Ansatz" (|12|). 
As shown in Appendix A and also on numerical examples in the main part of the paper this 
"Ansatz" guarantees rather good (quantitatively) approximation to an exact solution in the 
region of finite correlation lengths of fluctuations of short range order £, while in the limits 
of £ — > oo and £ — > our solution is exact. 

Calculation of the spectral densities shows that within both scenarios we can ob- 
tain rather attractive (in the sense of possible comparison with the experimental data on 
cuprates) picture of "destruction" of Fermi-liquid behavior on specific ( "hot" ) parts of the 
Fermi-surface, with persistence of Fermi-liquid state on the rest ("cold") part of the Fermi- 
surface. This non Fermi-liquid behavior is due to the strong scattering of electrons on 
fluctuations of short range order and, in general, is more visible with the growth of corre- 
lation length £. At the same time there are definite differences between these two scenarios 
which can be used, in principle, in the analysis of real experimental situation. In partic- 
ular, in the "hot spots" model (AFM-fluctuations) the pseudogap in the density of states 
is relatively weak (cf. Fig.9), while in the model of SC-fluctuations the pseudogap in the 
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density of states is more visible (cf. Fig. 13). The model of "dielectric" AFM-fTuctuations is 
more attractive even from simplest consideration of the phase diagram of Fig.l — pseudogap 
anomalies are mainly observed in the underdoped region and are apparently more intensive 
for systems which are closer to dielectric (AFM) state. It is obvious that precisely in this 
region we can expect more important role of fluctuations of "dielectric" (AFM) type and 
the growth of corresponding correlation length £. It is rather difficult to imagine why in this 
region of the phase diagram SC-fluctuations may become more important, this apparently 
must take place somewhere close to the optimal (corresponding to highest T c ) doping. Also 
in SC-scenario we have an obvious problem of characteristic scales (on temperature and 
energy) of pseudogap anomalies, which can not be solved within simple BCS-like theory, 
and requires some new microscopic approaches Our models are useful for the analysis 

within both scenarios of pseudogap formation irrespective of microscopic picture, because 
they are based on rather general (semi phenomenological) form of correlation function of 
fluctuations of short range order. 

Authors are grateful to O.V.Tchernyshyov for preliminary information on his analysis of 
one-dimensional model. 
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APPENDIX A: ANALYSIS OF ONE-DIMENSIONAL MODEL. 



Let us consider in more detail the use of the "Ansatz" QT2D to estimate the contributions 
of higher order diagrams. We shall limit ourselves to the analysis of one-dimensional model 
|15HT7| , because in one-dimension the problem is most serious |22|] . We are interested in the 
vicinity of Fermi "points" +pp ~Pf, while the Gaussian fluctuations of short range order 
scatter electrons by the momentum of the order of Q ~t 2pp from one end of the Fermi 
"line" to the opposite with scattering momentum values fixed with precision of the order of 
f^ 1 = k [0-[i7l| . We shall consider the linearized electronic spectrum: £- =t Vpp and 



to shorten notations put Vp — 1. Thus our system consists of "two types" of electrons - 
those moving to the "left" and to the "right". It is convenient to make our analysis in the 



coordinate representation p2| , when the equation of motion of electrons in our model takes 
the form 



("i-^^^fA-wT)*"' 1 ! < A1 > 

We limit ourselves with incommensurate fluctuations only when A*(x) 7^ A(x). Spinor 

^ = ^ ^ + ^ describes "right" and "left" electrons. Fluctuations A(x) are supposed to be 

Gaussian with < A(x) = > and < A*(x)A(x') >= \ A\ 2 exp(— k\x — x'\). Free propagator 
in frequency and coordinate representation has the form: 

G (ex) = i9(ea 3 x)sign(e)exp(iea 3 x) (A2) 

where 03 = +1 for the "right" particles, cr 3 = —1 for the "left". The particle passing the 
path of the length I produces a phase factor e l£l . During calculation of contribution of a 
given diagram it is convenient to change integration variables from coordinates of interaction 
vertices Xk to path lengths Ik, passed by the particle between scattering events |22fl . It is 



important to take into account the fact that these path lengths are not independent as for 
the given diagram the total particle displacement x — x' is fixed. The rules of the diagram 
technique to calculate G(e, x — x') are as follows : 



1. Electron line of length 1^ gives the factor —ie dk ^ £ ( ^ p \ 

2. Wavy (interaction) line connecting vertices m and n gives the factor: 
\A\ 2 exp(-K\x m - x n \) = \A\ 2 exp(-K\ E£=™(-1)%|)- 

3. Over all we must integrate from to 00. 

4. Integrate over p with a factor of e tp<yX ~ x '^ /2n . 

To calculate G{ep) just drop the last rule. From these rules we can see that the finite values 
of correlation length £ = k _1 lead to some damping of given transition amplitude with the 
displacement of the particle. The exact accounting of this effect is rather complicated but 
we can find some upper and lower bound estimates. Considering first the obvious inequality: 
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n— 1 n— 1 



«l E (-iW) > ezp(-« £ **) (A3) 



k=m k=m 



and using for the interaction line the r.h.s. expression in (|S3|), we overestimate transition 
amplitude damping (i.e. effectively overestimate k) for the given diagram. It is easy to see 
that the use of this approximation to calculate a given diagram for the Green's function leads 
(in momentum representation) to an extra ik appearing in every denominator of Green's 
function enveloped by an extra interaction line. This leads to an expression for the given 
higher order contribution of the form of dl2|) (cf. ||22|| ). For example, diagram shown in 
Fig. 14 yields (we assume e > 0, 5 = + ): 

AG(ep) = A 4 1 —— ( 1 —) l —— (A4) 

e — p + id^e + p + ine — p + 2ik e + p + ik j e — p + id 

which is analogous to (pl),(|Tl"l). On the other hand we can take the inequality: 

n—l n— 1 

exp(-K\ £ (-l) k k\) < exp(-K (-l) k - m h) (A5) 

k=m k—m 



and use for the interaction line the r.h.s. expression in flAq). This will lead to some underes- 
timation of damping effect in the given transition amplitude (i.e. effectively underestimate 
k). In particular, for diagram in Fig. 14 the contribution of interaction lines is: 



g— «J2g— fe— g— «J2g— Kpl— J2+J3) _ g— K(il+l3) (A6) 

In momentum representation this yields: 

AGM = A 4 ( — 1 —— 1 — ) 1 —— (A7) 

e — p + 10 x s + p + ik s — p + to e + p + ik / e — p + id 

Analysis of higher order diagrams shows that in this case contributions of all diagrams of the 
N-th order are equal and in the momentum representation we have ( "Ansatz of alternating 

G N {ep) = \A\ 2N 7 1 -^tt -, (A8) 



Then the whole series is easily summed analogously to the case of k = |13| , ^4]| and we 
obtain the Green's function in the form: 



5 It may seem that this choice for the interaction line contribution can even lead to the growth of 
transition amplitude in comparison with the case of k = and to appearance of some divergences, 
but this is not so. As we consider only incommensurate case here, where the interaction line 
envelopes only even number of interaction vertices (i.e. an odd number of l^), the choice of the 
sign in the exponent after removing the modulous is determined by dominance of even or odd 1^. 
This leads to effective damping of any diagram in higher orders. 
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G R (ep) = V N\G N (ep) = / dCe^- P . - s — ^ (A9) 

From this expression we can easily calculate the spectral density or the density of states: 
N(e) vpK 



N(E F ) ix J-J ip L d(e ~ (E* - e p ~ CW) 2 + (vFK) 2 (e ~ (A10) 

where we have restored vp explicitly. In Fig. 15 we compare densities of states for differ- 
ent values of k (correlation length), calculated using the "Ansatz of alternating k" and 
recurrence relations similar to (|l^) in one-dimensional model [15— 17| . We see that the re- 



sults are quantitatively close practically for all values of k. As we noted above our main 
"Ansatz" (|l^) , (|A4|) somehow overestimates the role of finite k, while the 'Ansatz of alter- 
nating k" ( |A7| ) underestimates it. Then it is guessed that the exact value of the density 
of states is in fact quite close to those obtained using these two types of approximation for 
the higher order diagrams. Analogous results can be obtained also for spectral densities. In 
fact this means that the results for the main physical quantities determined by one-electron 
Green's functions are not strongly dependent on the way the finite k enter the expressions 
for diagrams of higher order, but the main thing is to take account (at least approximately) 
of all diagrams of perturbation theory with different combinatorics. In principle this is not 
very surprising, as the main effect of pseudogap formation is due essentially to "backward" 
Q ~ 2pi?-scattering, which is accounted for exactly in the limit of £ — > oo, while the effect 
of finite k reduces to rather weak modulation of this random field, leading to the damping 
of its correlator and smearing of the pseudogap. 

Naturally, the "Ansatz of alternating k" can be written in the form of recurrence relation 
of the type of (|14}) also for two-dimensional, models discussed in the main part of this paper. 
For example in "hot spots" model we have: 

E k (e^ p ) = A 2 , , V{k) (All) 

l£n - ffc + ia k v k K - L fc+ i(£:„4p) 

where a k = 1 for odd k and a k = for even k. Other notations are given above in the 
main part. Data for the density of states obtained with the help of ( |A11| ) are shown above 
in Fig. 9 and confirm the our conclusions. Expression similar to ( |A11| ) can be easily written 
also for the model of SC-fluctuations. 

Let us stress that the "Ansatz of alternating k" is rather formal and is used only to 
show that this more or less arbitrary "approximation" , underestimating the role of finite k 
in diagrams of higher orders, leads to results which are quantitatively close those obtained 
with the "Ansatz" (|T2|),(A4), which overestimates this role. This last approximation used 



in Refs. ||T5|-|l^] and in the main part of this article has much deeper sense. As we have 
already stated above this approximation is exact in the vicinity of "hot spots" for the 
values of parameters of the "bare" spectrum t,t' and /i (topologies of the Fermi surface) 
which guarantee equal signs of velocity projections in "hot spots" connected by vector Q. 
Analogously in one-dimensional model it is possible to obtain the higher order contributions 



in the form similar to (|T2p or ( |A4p if we consider the model of correlator of fluctuations of 



short range order with the maximum at some arbitrary scattering vector Q which is much 
smaller than pp, so that (for large enough correlation lengths £) electrons are scattered by 
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fluctuations remaining always on one ( "left" or "right" ) branch of the spectrum. In this case 
expressions of the type of (K4) are exact. After that in final answers for diagrams of higher 
orders we can perform continuation to Q ~ 2pp of interest to us, as the only dependence on 
Q enters only via the "bare" electronic spectrum. Similarly we can achieve the same result 
changing appropriately the chemical potential fi (band filling). 
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Figure Captions. 



Fig.l. Schematic phase diagram of HTSC-cuprates [12|. For temperatures T < T cr 
there are well developed fluctuations of AFM short range order. For T* < T < T cr these 
fluctuations can be considered as static. 

Fig. 2. Model of the Fermi surface for HTSC-cuprates. Electronic states around the 
intersection points of the Fermi surface with magnetic Brillouin zone (shown by dashed 
lines) are strongly interacting with fluctuations of AFM short range order ("hot spots"). 

Fig. 3. First order self-energy diagram for the electron interacting with fluctuations of 
short-range order. 

Fig. 4. Second order sel-energy diagrams for the electron interacting with fluctuations of 
short range order. 

Fig. 5. Regions in spectrum parameters space where both "hot spots" exist (dashed 
region) and velocity projections in "hot spots " are of the same sign (doubly dashed region). 

Fig. 6. Fermi surfaces defined by the spectrum ([Sp, for different values of the chemical 
potential /i (band-filling) and parameter t' jt. 

(a) — case of t'/t = —0.6, when ([14]) is exact close to "hot spots": 

^/t=: (1)— -2.2; (2)— -1.8; (3)— -1.666. . . ; (4)— -1.63; (5)— -1.6; (6)— 0; (7)— 2, 
solution ( JT4] ) is exact in the vicinity of "hot spots" (velocity projections are of the same sign) 
for fi/t < —1.666 . . ., "hot spots" exist for \i/t < 0. 

(b) — case of t'/t = —0.4, characteristic of HTSC-cuprates, when ( |T4] ) is approximate: 
n/t =: (1)- -2.2; (2)- -2; (3)- -1.6; (4)- -1.3; (5)- 0; (6)- 2; (7)- 4, 

"hot spots" exist for -1.6 < fi/t < 0. 

Fig. 7. Energy dependencies of the spectral density in the "hot spot" (p x a/n = 
0.1666, p y a/iT = 0.8333) for different diagram combinatorics for the case of t'/t = —0.6, 
when fll4]) is exact: 

(a) — incommensurate case. 

(b) — commensurate case. 

(c) — combinatorics of spin-fermion model. 

Correlation length corresponds to the values of k: (1) — 0.01; (2) — 0.1; (3) — 0.5, 
A = 0.11 

At the inserts — energy dependencies of spectral density for different diagram combina- 
torics for na = 0.01 

(1) — at the "hot spot" p x a/n = 0.1666, p y a/n = 0.8333. 

(2) — close to the "hot spot" p x a/n = 0.1663, p y a/n = 0.8155. 

(3) — far from the "hot spot" p x a/ii = 0.0, p y a/n = 0.333. 

Fig. 8. Energy dependencies of the spectral density far from the "hot spot" (p x a/n = 
0.142, p y a/n = 0.857) for t'/t = -0.4, fj,/t = -1.3, which is approximately valid for HTSC- 
cuprates: 

(a) — incommensurate case. Dashed line — spectral density for incommensurate case ob- 



tained from ([All 
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(b) — commensurate case. 

(c) — combinatorics of spin-fermion model. 

Correlation length corresponds to the values of Ka: (1) — 0.01; (2) — 0.1; (3) — 0.5, 
A = O.lt. 

At the inserts — energy dependencies of spectral density for different diagram combina- 
torics for Ka = 0.01 

(1) — at the "hot spot" p x a/n = 0.142, p y a/n = 0.857. 

(2) — close to the "hot spot" p x a/n = 0.145, p y a/ix = 0.843. 

(3) — far from the "hot spot" p x a/n = p y a/7i = 0.375. 

Fig. 9. One-electron density of states for different diagram combinatorics ((a)-case of 
t'/t = -0.4, njt = -1.3; (b)-case of t'/t = -0.6,///* = -1-8): 

(1) — incommensurate case. 

(2) — commensurate case. 

(3) — combinatorics of spin-fermion model. 

(4) — in the absence of AFM fluctuations. 

Dashed line — incommensurate case, obtained from ( |A11|) . 
A/t = 1, correlation length corresponds to na = 0.1. 

At the inserts — one-electron density of states for commensurate combinatorics for: 
(1)— «a = 0.1; (2)— «a = 0.01 

Fig. 10. Self-energy diagrams in the model of SC-fluctuations: 

(a) — first order diagram with explanation of the meaning of interaction line (fluctuation 
propagator of Cooper pairs). 

(b) — second order diagram. 

Fig. 11. (a)-Energy dependence of the spectral density A(E,p) for the case of d-wave 
fluctuation pairing for different values of the polar angle 0, defining the direction of electronic 
momentum in the plane: 

(1)-0 = O; (2)— = tt/6. 

Correlation length corresponds to Vpn/A = 0.5 (full curve) and 0.1 (dashed). 
(b)-analogous dependence of the product f(E)A(E,p) (f(E)- Fermi function): 
(1)_ (j) = 0; (2)— = vr/6; (3)— = vr/4.83. 
Temperature (in Fermi function) T = 0.1 A, vfk/A = 0.5 . 

Fig. 12. Schematic picture of Fermi surface "destruction" by pseudogap due to fluctua- 
tional d-wave pairing. Dashed are the regions, where the spectral density is essentially non 
Fermi-liquid like. 

Fig. 13. One-electron density of states in the model of SC-fluctuations for the different 
values of parameter vpK/A: 

(a) — case of s-wave pairing. 

(b) — case of d-wave pairing. 

Curves are shown for the following values of Vpn/A: 
(1)— 0.1; (2)— 0.5; (3)— 1.0; (4)— 2.0. 

Fig. 14. Second order diagram for the correction to the Green's function in coordinate 
representation. 
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Fig. 15. One-electron density of states in one-dimensional model for the different values 
of Vpn/A: 

(1)— 0.1; (2)— 0.8; (3)— 1.2. 

Full curves — result of calculations using fll2D,fll4l) fTBH , dashed line — result of calcula- 
tions using ( |A10| ). 
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